Phase Separation in Mixtures of Repulsive Fermi Gases Driven by Mass Difference 
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We show that phase separation must occur in a mixture of fermions with repulsive interaction if 
their mass difference is sufficiently large. This phenomenon is highly dimension-dependent. Con- 
sequently, the density profiles of phase separated 3d mixtures are very different from those in Id. 
Noting that the ferromagnetic transition of a spin- 1/2 repulsive Fermi gas is the equal mass limit 
of the phase separation in mixtures, we show from the Bethe Ansatz solution that a ferromagnetic 
transition will take place in the scattering states when the repulsive interaction passes through 
resonance and becomes attractive. 



In the last few years, there have been considerable in- 
terests in strongly repulsive Fermi gases. Many of these 
studies were stimulated by the report of evidence of fer- 
romagnetism in the Fermi gas of 6 Li[l]. The possibility 
of itinerant ferromagnetism was first proposed by Stoner 
for electron gas (2]. The idea is that if Coulomb repulsion 
increases faster than kinetic energy as a function of in- 
creasing density as indicated by Hartree-Fock estimates, 
the system will turn ferromagnetic and form ferromag- 
netic domains at sufficiently high densities. However, 
Hartree-Fock approximation overestimates repulsion en- 
ergy. So far, itinerant ferromagnetism has not been found 
in metals. 

Itinerant ferromagnetism had also been predicted for 
strongly repulsive Fermi gas based on perturbative and 
mean field calculations 0E] prior to the MIT experiment 
PQ . However, such approaches are known to be unreliable 
in strongly interacting regime. In fact, no ferromagnetic 
domain were observed in ref.[T], and newer experiments 
have established that the repulsive Fermi gas of 6 Li is 
non-magnetic [5] . This shows that within the range of in- 
teraction studied experimentally, the repulsion in Fermi 
gases does not dominate over the kinetic energy as mean 
field calculations show. Yet Stoner's idea should work. It 
should also be applicable to systems such as binary mix- 
tures of Fermi gases, (referred to as Fermi- Fermi mixture 
for short). In that case, the analog of ferromagnetic tran- 
sition of spin- 1/2 Fermi gas is phase separation of the two 
fermion species. 

Phase separation of Fermi-Fermi mixtures has been 
studied in ref . [6] using mean field approximation and per- 
turbation methods. It is found that a 6 Li- 40 K mixture 
will phase separate in the strongly interacting regime. 
Since the same mean field theory predicts a ferromag- 
netic transition in a spin- 1/2 Fermi gas in the strongly 
interacting regime £Q which is not observed [5], it raises 
the questions about whether increasing repulsion can in 
fact cause a Fermi-Fermi mixture to phase separate as 
mean field theory predicts. 

In this paper, we would like to point out that phase 
separation in a Fermi-Fermi mixture can always be in- 
duced by increasing the mass ratio of the two fermion 



species, but not necessarily by increasing repulsion. The 
reason is that the kinetic energy cost for phase separa- 
tion can always be reduced to zero by increasing the mass 
ratio, thereby falling below the repulsion energy, render- 
ing the Stoner argument valid. On the other hand, since 
the density regime for strong interaction is dimension de- 
pendent, the phenomena of phase diagram changes sig- 
nificantly with dimensionality. Since the ferromagnetic 
transition in spin-1/2 systems (which leads to formation 
of domains) is the equal mass limit of the phase sepa- 
ration of Fermi mixtures, it is useful to unify these two 
phenomena with a global phase diagram as a function of 
mass ratio and interaction. In the Id case, we shall also 
show from exact result that an "upper-branch" spin-1/2 
Fermi gas will turn ferromagnetic as it passes through 
resonance to the side of attractive interaction. 

(A). A theorem on mass-difference driven phase sep- 
aration: The energy density £hm of the ground state of 
a homogenous mixture of heavy and light fermions with 
masses (mi,, mji) and densities (nj,, nn) is 

£hm = £l + £h + £lG [ — - , nl d a, — ] , (1) 

where £ L (H){n L (H)) = A dn ( ^^ /d /m L{H ) is the energy 
density of the ideal gas of the light (heavy) fermions, d is 
the dimensionality, and Ad is a constant. The last term 
U = £t,G is the interaction energy in units of £j,, and 
G is a dimensionless function of the variables displayed, 
"a" is the length scale associated with the interaction. 
In 3d, a is the s-wave scattering length a s in the pseudo- 
potential 
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We set h = 1 in this paper. By applying harmonic con- 
finement along the axial (with frequency lv z ) or the trans- 
verse (ui) direction, the system can be reduced to a 
quasi 2d or a quasi Id system. For quasi 2d systems, 
a is related to the binding energy as Cb = l/(2ma 2 ), 
where = ^u z e^ a ~l a s ; Gz = yj 1/ (muj z ) is the con- 
finement length and B w 0.915J7 . For quasi Id sys- 
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tems, a — — — C) where a± = y/l/(muj±) and adiabatic theorem, dEh m /d( > 0, where £ ar e — l/a[H], 

C ~ 1.46|8 . In all dimensions, the energy satisfies the In a [10, and a[TT] respectively for 3d, 2d and Id systems. 




FIG. 1. Figure 1A, IB, and 1C are the phase diagrams for a 3d, 2d, and Id Fermi-Fermi mixture with Nl = Nh in a volume 
V. kp = (67r 2 n) 1//3 and -nn respectively for 3d and Id, with n = Nl/V = Nh/V. To the right of the vertical blue line, the 
mean-field interaction energy is less than half of total kinetic energy for a homogenous mixture (L,H), and the system is weakly 
interacting deeper in that region. The gray dashed-dot lines indicate the case of a Li-K mixture with niL/ran = 6/40. Figure 
ID, IE, and IF are the phase diagrams of a 3d, 2d, and Id Li-K mixture in chemical potential plane for weak interactions. 
jiH, flL are scaled by (fn d g 2 ) 1 ^ 2 ~ d ' > in 3d and Id, and 1/m in 2d. The red dashed-dot lines in ID and IF represent trajectories 
for the density profiles of a trapped system, corresponding to (a-d) in Fig. 2, with the squares denoting the chemical potentials 
at the trap center. From Figure 1A to IF, the black (orange) solid lines represent the 1st (2nd)-order boundaries with (without) 
density discontinuity. In Figure IB, the boundary is given by the function j c ( m ti™fl) — 2-7T / \fmtmH ■ In IE, the two solid 
orange lines are the boundaries for interaction g < g c , with two slopes (g/g c )y/m,H /m^ and (g c /g) \J tub /ttil respectively. 
When g > g c , the two boundaries merge into one (shown by dashed line) with slope \ frniTfmL,- 



To find the condition for phase separation, we consider 
a system with Nl and Nh fermions in a volume V . Let 
us define 

m L /m H = x, N H /N L =~{, (3) 
the total energy of the homogenous mixture is 



E h 



y£ L (n L )(l+7 Q x + G), a = l + 2/d. (4) 



Next, we consider the fully phase separated state. Let Vh 
and Vl be the volumes of the heavy and light fermions, 
Vh + Vl = V. The ratio Vh/Vl is determined by equat- 
ing the pressure P of these two separated gases. Since 
the pressure of an ideal gas is proportional to its energy 
density, P — 28 /d, we have £l(ti' l ) = £y(n^), where 



N H (L)/V H (L)- This gives Vh/Vl 



l/a 



The total energy of the phase separated state is Eps 
V H £ H {n' H ) + V L S L (n' L ) = V£ L (n L )(V/V L )^ +d ^ d , or 



E 



PS 



V£ 



(5) 



The phase separated state will have lower energy if Eh m - 
Eps > 0, or 



I(x) = G{x) - \(l + jx^ a ) a - 1 - 7 ° 



> 0. 



(6) 



When the mass ratio is sufficiently small such that 
7 1 / Q a; -C 1, hence x a < x, Eq.|6| becomes 

(7) 



I(x) = G(0) - a-fx^ a + 0(x,x 2 /°) > 0, 



where G(0) > is the repulsive interaction energy in the 
limit when m# — > oo|12). Eq.Q can always be satisfied 
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for sufficiently small x, thus proving that phase separa- 
tion must occur for sufficiently large mass difference. 

Corollary: Because of the adiabatic theorem, if a 
mixture with mass ratio m^/m^ phase separates at a 
given interaction £, it will continue to phase separate at 
stronger interactions, i.e. at a larger £. 

(B). Phase diagram: To be concrete, we shall consider 
the phase diagram for weak repulsion. In this case, the 
interaction energy U in Eq.Q is much less than the ki- 
netic energy, and is accurately described by the mean 
field form U = gn L n H , where g = ^ in 3d, in 
quasi 2d, and in quasi Id. The energy density £hm, 

chemical potential (/if,, (J-h) for light and heavy particles, 
and the pressure of the homogenous mixture are 

£hm(n L ,n H ) = £ L (n L ) + £n(nn) + gn L n H , (8) 

, v d£ L{H) {n L{H) ) 
^L(H){n L ,n H ) = \-gn H ( L ), (9) 

P(n L ,n H ) = pt L n L + ^ H n H - £{n L ,n H ). (10) 

For a system with Nl light fermions and Nh heavy 
fermions in a volume V, the possible configurations 
are: (a) fully phase separated state (PS), denoted as 
(LlkH); (b) coexistence of a homogenous mixture and 
a single phase, denoted as (L,H)kL or (L,H)SzH; 
(c) coexistence of two homogeneous mixtures with dif- 
ferent densities (n' L ,n' H ) and (n'l,n H ), denoted as 
(L' , H')&z(L" , H"); and (d) a single homogenous mix- 
ture (L,H). To determine the presence of these 
phases, it is sufficient to consider the general situation 
(L 1 , H')Sz(L", H"), since it covers all other cases. For 
example, the state (LSzH) corresponds to n' H = n" L = 0. 
The state (L,H)k,L corresponds to n' H — 0, and the 
state (L, H) corresponds to n'[ = n H = 0. 

Let us denote the particle numbers in the two coexis- 
tent homogenous mixtures as (N' L ,N' H ) and (N' L \N H ), 
and their volumes as V and V" , where N' L + N' L ' = Nl, 
N' H + N H = Nh ; and V + V" = V. Minimizing the to- 
tal energy with respect to N' L , N' H and V will determine 
the type of configuration ((a) to (d)) for given interac- 
tion, from which one can construct the phase diagram. 
In the cases of (b) and (c), the minimization procedure 
leads to the Gibbs phase rules, i.e. two coexisting phases 
have the same pressure and the same chemical potential 
for the same atomic species. 

The phase diagram of a 3d mixture as a function of 
chemical potential and interaction has been studied in 
ref.[B]. Here, we examine the phase diagram in different 
dimensions, focusing on the effect of mass ratio. Figure 
I A, IB, and I C show the phase diagrams for a 3d, 2d, and 
Id mixture with Nl — Nh in a volume V. In all cases, 
the system is fully phase separated (L&H) for sufficiently 
small mass ratio mi/m ff . For both Id and 3d, there is 
a range of mass ratio (for given interaction) in which the 
system consists of two different phases in equilibrium, 



((L,H)kL for Id and (V ,H')k(L" ,H") for 3d). This 
feature is absent in 2d [13]. 

Note that the phase boundaries in Figure 1A to 1C are 
derived from the weak coupling form of the interaction, 
they are therefore inaccurate in the strongly interacting 
region. However, the corollary in section (A) guarantees 
that the system will phase separate in the strongly inter- 
acting region over a range of mass ratio wider than that 
in the weakly interacting region. 

(C) Ferromagnetic transition of Id spin-1/2 Fermi 
gas at resonance: The phase diagram for Id Fermi- 
Fermi mixture is not only constraint by the results in 
the weakly interacting regime, but also by the exact 
Bethe Ansatz solution along the line mi/mjy = 1[14). 
which is a spin-1/2 repulsive Fermi gas with interaction 
g J2i>j $( x i ~ x j)i were g — — (ma) -1 . Because of the in- 
tegrability of this system, there are two classes of eigen- 
states: one where all quasi-momenta are real, i.e., all 
particles are in scattering states, (denoted as class (i)), 
and one that contains at least one pair complex conjugate 
quasi-momenta, i.e. with at least one fermion bound pair, 
(denoted as class (ii)). Repulsive Fermi gas is in class (i). 
It is stable without decay into class (ii)|T5J. System in 
class (i) will be referred to as in the "upper branch" . 

Experimentally, one can tune the interaction from the 
repulsive regime (—1/g < 0) to the strongly repulsive 
regime near resonance (—1/g = - ), and then to strongly 
attractive regime (—1/g = + ). The ground state of a 
repulsive (—1/g < ) spin-1/2 Fermi gas with equal spin 
population is a spin-singlet according to the Lieb-Mattis 
theorem|16j. Its energy is a monotonically increasing 
function of —1/g according to the adiabatic theorem. At 
resonance, —1/g — 0, the spatial wavefunction of the 
ground state is identical to that of a fully spin polarized 
Fermi gas up to a sign (which changes in various regions 
in configuration space). As a result, its energy E(0) is 
given by that of a fully spin polarized state with huge spin 
degeneracy [TJ] - all spin configurations including the spin 
configurations (a) to (c) mentioned above are degenerate, 
with H and L now labeling the two spin species. This 
means that the two phase boundaries in Fig.lC will con- 
verge to the equal mass point tol/toh = 1 at resonance. 
Crossing the resonance to the attractive side, —1/g > 0, 
the energies of all spin states continue to increase, hence 
E( — l/g > 0) > E(0); except for the largest spin state 
which remains at E(0) regardless of interaction. The sys- 
tem will therefore make transition to this maximum spin 
state. In the presence of a small magnetic field with a 
weak field gradient that destroys spin conservations, the 
spin in the ground state will align along the field direc- 
tion. 

The equal mass line in Fig.lC shows that the sys- 
tem will not phase separate (or turn ferromagnetic) even 
when repulsion is made the strongest (i.e. at resonance). 
Similar situation occurs in 3d. There, experiments show 
that the spin-1/2 Fermi gas is not ferromagnetic (i.e. 
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no phase separation) at all strongly interacting regime 
examined [5]. Instead, by increasing the mass difference, 
the repulsive mixture can definitely phase separate in any 
dimension, as guaranteed by the theorem in Section (A) . 
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FIG. 2. Density profiles of a trapped Li(Light)-K(Heavy) 
mixture in 3d((a) and (b)) and Id((c) and (d)), corresponding 
to the trajectories (a) to (d) in Fig. ID and IF. The densities 



(nH,riL) have been normalized by 



the density of light 



atoms (Li) at the trap center for non-interacting system in the 

same setup. The position r is scaled by a; lo = ^I/^lujl), 
the confinement length of light atoms, (a) and (b) are with 
the same particle numbers (N H ,N L ) = (10 5 )(1.47, 6.9) and 
the same trapping frequency ratio ujh/^l = 0.3, but with 
different interaction strengths a/a^o = 0.055(a), 0.065(6). 
(c) and (d) are with the same (N H ,N L ) = (10 4 )(1.45, 0.51) 
and the same interaction — a;, D /a = 157T, but with different 
Wfl/wi = 0.5(c), 0.2(d). 



(D) The density profile in a trap: Phase separation can 
be easily identified in the density profile of a trapped gas, 
which can deduced from the phase diagram in the chem- 
ical potential /i^ — fin plane. To illustrate the possible 
classes of phenomena, we again consider weak interac- 
tions. 

Given (i^ and three phases are allowed: the sin- 
gle component Fermi gas (L), (H), and the homogenous 
mixture (L,H). To express the pressure of a homogenous 
mixture Phm as a function of fir, , fin , we invert Eq. Q to 
obtain nj, and nu as a function of \xl and and then 



substitute them into Eq.( 10 ). The pressure of (L) or (H) 

is Pl(h)((i>l, Mh) = B d' n H(H)^ 1 L{H) '> wnere B d is a con- 
stant. The phase boundary for the full phase separation 
is Pl(hl) = Ph(Hh), or 



H H /li L =/3, = (m L /m H ) d ^ d+2 \ (11) 



The phase boundary between the mixture {L, H) and 
L (or (H)) is obtained by equating Phm{^L,^H) — 
Pl(h)(hl(h))- The phase boundaries for the 3d, 2d, and 
Id mixtures are shown in Figure ID, IE, and IF respec- 
tively. Within the region of homogenous mixture, the 
inversion of Eq.(|9]) may yield several solutions of densi- 
ties (say, (n' L ,n' H ), (n'^n'jj)) for given chemical poten- 
tials (/iL,/i//), which in general have different pressures. 
The thermodynamic state is given by one with highest 
pressure. In the 3d case, the homogeneous mixture is 
contained within the "bubble" in Figure ID. Within this 
region, the thermodynamic state is unique except on the 



line that is an extension of the boundary Eq.(ll) where 
two states (with densities (n' L ,n' H ), (n'^n'jj)) have iden- 
tical chemical potential and pressure. The difference in 
density across the boundary can be obtained by noting 



that the thermodynamic relations (Eq.([8|) to (10)) of the 
mixture is invariant under the transformation 



i' L = (3n Hl n H -> n' H = f3 1 n L , (12) 



The density discontinuities across this line are then 
Aul — ftnjj — Ul, Anjj = j3~ 1 nr / — nu, which implies 



An L /An H = -p. 



(13) 



In Fig. 2a to 2d , we show the density profiles of the 
3d and Id mixtures in a trap constructed by applying 
local density approximation (LDA) to the equation of 
state n L ( H) (v) = n L ^ H) (fi L - V L (r),fi H - V H (r)), where 
Vj,(r) and Vff(r) are the potentials experienced by the 
light and heavy particles. Moving from the center of the 
trap to the surface of the cloud corresponds to following 
the trajectories indicated in Fig. ID and IF. Fig. 2a and 
2b show the density profiles of a 3d mixture at different 
interaction strengths. The discontinuities of the densities 
follow Eq.(13). Fig. 2c and 2d show a Id mixture under 



different trapping potentials. 

Two features of the density profiles should be empha- 
sized. Firstly, the density profiles of a 3d mixture differ 
significantly from that of the Id mixture, (see Fig. ID 
and IF). Phase separation takes place in the outer part 
of the atom cloud in Id but in the inner part in 3d. 
This is because the strongly interacting regime occurs 
in the low (high) density region in Id (3d). Secondly, 
in Fig. 2a-2d, we note that n^H) can increase with 
r. This is different from the single component case, 
where dn/dr < 0, due to the fact that dn/dfi > 
as demanded by thermodynamic stability. In the mix- 
ture case, stability against density fluctuation requires 



Det(M) > 0, where M i:j = 
and H. We then have dni/dr 
Ah -9 
-9 A L 



dfii/drij, and i,j = L 
= (M~ 1 ) i jdfj,j /dr, where 



M 



Det _1 (M) 



.4 



L(H) 



dfi L 



> 0. 



dn L(H) 

That dn L ( H j/dr can be positive or negative is because it 
is made up of two terms. If driL /dr > 0, it is easily shown 
from stability condition (A^Ah > g 2 ) that dnn/dr < 0. 
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Thus one can have at most one species with a positive 
density derivative. 

Finally, we point out that the mass-difference mecha- 
nism for phase separation discussed in Section (A) will 
also be effective for mixtures with more atomic species; 
and the original Stoner idea will be illustrated most ready 
in the large mass difference limit. The features mentioned 
in this section, and the ferromagnetic transition of an 
upper branch spin-1/2 Fermi gas as implied by Bcthe 
Ansatz should all be experimentally observable. 
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